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Abstract 

It has been predicted that dilepton production may be used as a quark-gluon plasma 
probe. We calculate the rapidity distributions of thermal dileptons produced by an 
evolving quark-gluon plasma assuming a longitudinal scaling expansion with initial 
conditions locally determined from the hadronic rapidity density. These distribu- 
tions are compared with Drell-Yan production and semileptonic charm decays at 
invariant mass M = 2, 4, and 6 GeV. 
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One of the challenges in ultrarelativistic heavy ion physics is to differentiate the ther- 
mal signatures of a quark-gluon plasma from particles produced in the initial nucleon- 
nucleon interactions. In order for dilepton production to be a useful plasma probe one 
must then distinguish between thermal production and the significant 'background' dilep- 
ton production from initially produced Drell-Yan and charm pairs. We focus on the 
rapidity distributions of lepton pairs in the mass region between 2 < M < 6 GeV where 
thermal production could be competitive with initial dilepton production. We also show 
that thermal charm production is significant. 

The rapidity distribution of Drell-Yan pairs of invariant mass M in central (impact 
parameter 6 = 0) AA interactions, averaged over protons and neutrons, is defined as 
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where K = 2 accounts for higher-order corrections, v(x) = u v (x) + d v (x) is the valence 
quark distribution, and s(x) = 2(u(x)+d(x)+s(x))/6 is the average sea quark distribution 
in a nucleon assuming a flavor symmetric sea. We have used the new lowest-order set 
of structure functions, Duke-Owens 1.1 in our calculations. The fractions of the 
initial projectile and target momenta carried by the interacting partons are given by 
x a = Me y /yfs and x b = Me~ y /^/s, respectively. 

The rapidity distributions from Drell-Yan production in a central Au+Au collision 
are shown in Fig. 1(a) for RHIC energies and 1(b) for the LHC. At RHIC, the rate rises 
slightly for M = 2 GeV and stays fairly constant over three units of rapidity even for 
M = 6 GeV. It drops off steeply at the kinematic limit. At the LHC, the rate rises 
somewhat over four units of rapidity for all masses. This effect is most prominent at the 
LHC for M = 2 GeV where x b < 0.00032 and x a > 0.00032. The sea-sea term in eq. (1) is 
the dominant contribution and increases slightly since x b remains near the small x limit 
of the structure functions while x a increases with y. The first sea-valence contribution in 
eq. (1) increases with y but the second decreases because a small decrease in x b leads to a 



fast decrease of x b v(x b ). Thus the total Drell-Yan rate increases with y by ~ 30 percent. 
The charmed hadron production cross section in a nucleon-nucleon collision is [§, || 
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where subscripts a and b represent the incoming partons, 1 and 2, the charmed quarks, 
and 3 and 4, the charmed hadrons. The function D(z) describes the fragmentation of 
the charmed quarks into hadrons. In ref. it was shown that charm hadroproduction 
is best described by the assumption that the charmed quark loses no momentum dur- 
ing hadronization, i.e. D(z) = 5(1 — z). Given this, e.g. p^ = pu 3 , but y 1 7^ y 3 . (A 
fragmentation function assuming some momentum loss results in a more central rapidity 
distribution.) The convolution function H(x a ,x b ) is 

H(x a ,x h ) = ^2(x a q nf (x a )x b q nf (x b ) + x a q nf (x a )x b q nf (x b ))^\ qq - + x a g(x a )x b g(x b )^\ gg , (3) 

where q(x) and g(x) are the parton distributions |J. For the explicit subprocess cross 
sections, see e.g. ref. [§J]. 

We assume no intrinsic Pt for the light quarks and use four-momentum conservation 
to integrate over x a and x b . Then the number of charmed hadrons produced in a central 
AA collision is 
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where M D -^ and y D p are the mass and rapidity of the DD pair. The uncertainty in the 
magnitude of the charm production cross section from a lowest-order calculation results 
in K ~ 2 — 3. The shape of the rapidity distribution is relatively unaffected by next-to- 
leading order corrections ||. Since pA measurements || [!], || indicate a ~ 0.9 — 1, we 
use a = 1 to give an upper bound on initial charm production. Using eq. (4), we can 
determine the rapidity distribution of lepton pairs from the decay DD — > l + l~ + anything 



and from the average branching ratio B(D/D* — > l + + anything) ~ 12%. (B(D° — > 
/+ + anything) = 7.7% and £(£>+ -> /+ + anything) = 19.3% while D°* -> D°X and 
Z) + * — ► D + /D°X where X is either a photon or pion.) Given a DD pair with mass and 
rapidity distributed according to eq. (4), the decay to dileptons is performed. The pair 
is first assigned a transverse momentum and polar angle using a parametrization of the 
charm pair correlation data of ref. ||. Then the two D's are semileptonically decayed 
according to the lepton energy spectrum reported in ref. [IU|. At high energies and for 



large nuclei more than one DD pair may be produced so that the detected charm decay 
leptons may be uncorrelated. 

The charmed hadron rapidity distributions are shown in Figs. 2(a) and (b) for Au+Au 
collisions at RHIC and LHC (solid curves) along with the resulting lepton pair distribu- 
tions at M = 2 (dashed), 4 (dot dashed), and 6 (dotted) GeV (scaled up by 100 to fit on 
the same plot). Charm decays are copious at these energies (cr D -p ~ 190 fib at RHIC and 
1.95 mb at LHC using m c = 1.5 GeV and K = 2, compatible with the next-to- leading 
order results ||). Measurements suggest that charm decays are important in the dilepton 
spectrum at M = 2 GeV at SPS energies JTTJ, O. Fermilab experiment E789 at 800 GeV 



observed that the dilepton continuum at M < Mj/^ is much greater than expected from 
Drell-Yan alone, indicating that the dominant background may be due to charm decays 



i~3|j . A judicious choice of kinematic cuts can reduce the charm acceptance relative to 
that of Drell-Yan. Drell-Yan pairs are produced with equal but opposite transverse mo- 
mentum with pt < M while the transverse momenta of the leptons from charm decays 
are only weakly correlated. This is discussed in more detail below. 

We have not included shadowing of the nuclear structure functions in our calculation. 
Shadowing would lead to a depletion of both the charm and Drell-Yan cross sections 
in the very low x^ region. The effect is strongest at small y since both the beam and 
target distributions are shadowed, leading to a faster increase of the Drell-Yan rate with 
y. Aside from the uncertainties in the rates due to shadowing, it is important to note that 
the nucleon structure functions have not yet been measured in the LHC x region and are 



also not well known for RHIC. HERA will study the region 1CT 5 < x and Q 2 > 10 GeV 2 . 
Since extrapolations of the measured structure functions to the low x region differ widely, 
the HERA results could, by favoring some models over others, lead to significant changes 



in the calculations of ctdy and cr D p. For example, using KMRS BO Hl4[ , a higher-order 
set of structure functions with a larger sea quark contribution at low x, increases the 
Drell-Yan cross section by a factor of two. Early results from HERA favor sea quark and 
gluon distributions that increase as x — ► 0, suggesting that, particularly at LHC energies, 
QCD calculations are quite sensitive to the choice of structure functions ]T5[ . 

Thus far we have been considering the expected background to the thermal dileptons, 
now we discuss the signal itself. At RHIC and LHC it is reasonable to expect the longitu- 
dinal expansion to follow the scaling law, v z — z/t, so that the initial density is related to 
the final multiplicity distribution by n^i = (dN ' j ' dr\) j '(irR 2 ^. The temperature is therefore 
a function of 77, T 4 = Tiiji). The multiplicity distribution is parametrized as 

dN (dN\ . 2/0 2 , , CN 
-=^-j o exp(V/2. 2 ), (5) 

where (dN/dr]) is the total multiplicity at rj = and a is the width. At RHIC energies, 
< 6 and we choose (dN/dr]) Q = 2000 while for the LHC \rj\ < 8 and we take (dN/di]) = 
5000. Assuming that the Gaussian shape extends to the kinematical limit, the energy sum 
rule is saturated for o ~ 3 at RHIC and ~ 5 at LHC. This form as well as the scaling 
properties of the longitudinal velocity must, of course, break down at rapidities close to 



the phase space limit [[16] where the densities become too small for thermal production 
to be observed. 



The dilepton emission rate is [17, 18 



dN a 2 
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where M, p?, vn^i and y are the mass, transverse momentum, transverse mass, and 
rapidity of the pair and rj is the rapidity of the medium with temperature T. In a plasma, 
Fq — e q) while in an equilibrium hadron gas, assuming tttt — > p — >• is the dominant 



channel, Fjj = j^TTip/ (( m p — M 2 ) 2 + m 2 T 2 ). The mixed phase has a fractional contribution 
from each. After integration over px, t, and transverse area, the rapidity distribution for 
a given dilepton mass is 
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where x = M cosh(i] — y) /T, Sj is the initial entropy density and sq and s# are the entropy 
densities of the plasma and the hadron gas at the transition temperature T c . 

We neglect the transverse expansion and describe the properties of the matter in the 
longitudinal direction only. Transverse expansion mainly affects the hadron gas phase 



and the later part of the mixed phase [[1^, |20[. If the transverse expansion is included 



the contribution from the hadron gas phase is negligible for M > 2 GeV except at large 
rapidities where the initial density is too low to produce a plasma. 

As defaults, we assume a three-flavor plasma and a hadron gas of massless pions with 
T c = 200 MeV and T dec = 140 MeV. Then 

S . = S ^S 7 ? Si > S Q or Si < Sh 

\ foSQ + (1 - f )s H S H < Si< Sq 

where 7& is the number of degrees of freedom with jh = 3 and 7q = 16 + 21n//2 and 
rii = Si/3.6. The beginning of the mixed phase occurs at 



1 fm Ti = T r . 



The initial plasma content of the mixed phase is 



1 T > T c 
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The beginning of the hadron phase is then 
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where r = sq/ Sh = Jq/jh- There is no thermal contribution if Tj < T dec . 

We have used three different hypotheses to fix T, and Tj. First we assume a fixed 
initial time, Tj = 1 fm. We then define dN/dr]\Q = nR\nQTi, dN/drjln = 7cR\nHTi, and 
dN/drj\ dec = iiR\ndec T i- If dN/drj > dN/drj\Q the system begins in the plasma phase and 

UN T 3 1 1/3 
[ drj nR z A n Q Ti J 

If dN/dr]\H < dN/dr] < dN/drj\q, the system begins in the mixed phase. If dN/drj\& cc < 
dN/drj < dN/drjln, it begins in the hadron gas with 

UN T 3 1 1/3 
Tl Th \ dr] nR 2 A n H Ti J ' ^ 

There is no contribution if dN/drj < dN/drj\a e c- The maximum initial temperature, Tj )max , 
obtained with a fixed initial time is the lowest of all three hypotheses. In this case we 
find T iimax ~ 255 MeV at RHIC and 350 MeV at the LHC. 

As an alternative way of determining the initial time Tj and temperature Tj, we use 
the bound Ti{Ej) ~ 3rjTj = C > 1 from the uncertainty relation. We consider two cases, 
C = 3 and C — 1, where the latter gives an absolute lower bound on Tj and an upper 



bound on Tj. Resent results from the parton cascade model 0] and estimates of gluon 
thermalization |22| give some support for a short formation time. In this case, the initial 
temperature of the plasma is 

T i dN 3r ° r do) 

If Tj < T c , then the formation time is again fixed to be Tj = 1 fm and the system begins 
in the mixed phase or hadron gas, according to the multiplicity, as before. When C = 3, 
we find Tj, max ~ 300 MeV at RHIC and 465 MeV at the LHC. Using C = 1, we have 



Tj imax ~ 515 and 810 MeV respectively, similar to the results of ref. p3 |. 

In Fig. 3 we show thermal dilepton production at RHIC for M = 2 GeV (a) and (b), 
M = 4 GeV (c) and (d), and M = 6 GeV (e) and (f). The full evolution (plasma, mixed 
phase, and hadron gas) is shown in (a), (c), and (e) while the plasma contribution alone is 
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given in (b), (d), and (f). The thermal dilepton distributions arising from the fixed initial 
time with Tj from eq. (8) (solid) and the uncertainty bound from eq. (10) with C = 3 
(dashed) and C — 1 (dot-dashed) are given. Dilepton production by the plasma is most 
important at central rapidities where the initial temperature is greatest. Including the 
mixed phase and the hadron gas broadens the rapidity distributions. Without transverse 
expansion, the lifetime of the hadron gas phase at central rapidities, where the initial 
density is high, becomes very long (of order 80 fm when 7] ~ Tj max ) and its contribution 
to the thermal pair rate dN/dMdy is thus unphysically large. At higher rapidities, the 
calculated hadron contribution is more physical since the initial density is too low to form 
a plasma and the effects of transverse expansion remain small. In Fig. 3(a), (c), and (e) 
the thermal dilepton distributions merge at higher rapidities since 7* = 1 fm when Tj < T c . 
Later time contributions are less important when C — 1 because of the high T i max . The 
distributions from plasma alone, Fig. 3(b), (d), and (f), are narrower and have a reduced 
yield. Since it is difficult to produce high mass pairs when Tj < T c , at M = 6 GeV the 
yield at y = is unchanged in Fig. 3(e) and (f). The results for the LHC are given in Fig. 
4. Due to the increased T i max , the contribution from the plasma dominates the thermal 
emission for M > 2 GeV at this energy. To show the most optimistic results of thermal 
dilepton production, we will compare the C — 1 curves with the other contributions to 
the dilepton spectrum. One must keep in mind, however, that the other choices of initial 
conditions result in a much smaller thermal yield for M > 2 GeV. 

We have also varied our inputs to see how the rapidity distributions are affected. 
Increasing the hadron degrees of freedom reduces the mixed phase and hadron gas contri- 
butions since their lifetimes are decreased. Changing to a two-flavor plasma increases the 
initial temperature considerably. For 7* = 1, T iyUiax grows to 285 MeV at RHIC and 386 
MeV at the LHC. When C = 3, we have 340 MeV and 540 MeV while when C = 1, we 
find 590 and 930 MeV respectively. Setting T c = 140 MeV and T dcc = 100 MeV increases 
r m but results in a decrease of the mixed phase and hadron contributions due to the larger 
x c and Xdec (lower T c and T dec ) in eq. (7). 
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If C = 1 represents the more physical situation, one might expect a significant amount 
of thermal charm production since Tj jmax ~ (1/3 — l/2)m c . The importance of thermal 
charm production has previously been suggested |22| , For completeness, we also esti- 
mate the rapidity distribution of dileptons from thermal charm production. The thermal 
cc production rate is 

where v re \ is the relative velocity of the thermal partons, v re \ = J {k a ■ k^) 2 — m^/(E a E b ). 
We use a Boltzmann distribution for the incoming partons, as before. Following ref . |24|] , 
we write the charm production cross section as 

a c - c (M) = 7,<t^«(M) + lg a gg ^{M) , (12) 

where the quark and gluon degeneracy factors are 7 g = 3 x (2 x 3) 2 for three quark flavors 

and 7 9 = (2 x 8) 2 /2 (the factor of 1/2 is to prevent double counting) p4j . The cross 

sections are given in ref. [PJ. We consider only production from the plasma phase for 

this case since the later stages of evolution do not have sufficiently high temperatures to 

provide significant thermal charm production. We have not included charm fragmentation 

since this has only a minor effect on the rapidity distribution as compared to that of the 

source. Therefore, we may simply write 

dN 3a c - c (M) 2 , [Y „. 2 x 4 + 5z 3 + 15x 2 + 30:r + 30 

>i$0(si - s Q ) / d^Tff -g- e X: ■ (13) 
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The resulting distributions are similar in shape to the thermal dileptons but broader 
than the initial charm production. Thermal charm production is comparable to that in 
the initial nucleon-nucleon interactions. In Figs. 2(c) and (d) we compare the thermal 
and initial charm production. The mass distribution is much steeper for the thermal 
production, resulting in a smaller decay contribution at higher masses. For comparison, 
at RHIC, the M = 2 GeV curve is multiplied by a factor of 1000 while the M = 4 and 6 
GeV results are multiplied by 10 4 and 10 5 respectively. At the LHC, the multipliers are 
100 for M = 2 GeV, 1000 for M = 4, and 10 4 for M = 6 GeV. 
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Our combined results for lepton pairs with M = 2, 4, and 6 GeV are given in Figs. 5 
and 6 at RHIC and LHC respectively. The Drell-Yan and initial charm contributions are 
given by the dashed and dot-dashed curves. The thermal dileptons arising from direct 
production over the full time evolution and charm decays produced by plasma only with 
C = 1 are given by the solid and dotted curves. The thermal dileptons have the broadest 
distribution of all the contributions at RHIC while the decays of initially produced charm 
have the narrowest. The slopes are also somewhat different near the central region. At 
M = 2 GeV, the thermal contributions are above the Drell-Yan. At the LHC, the thermal 
rate remains compatible with the Drell-Yan yield and all the slopes are similar. However, 
at y > 3 the thermal and charm decay contributions are decreasing while the Drell-Yan 
rate is still increasing. 

Although the initial charm contribution is very large, it may be possible to separate 
it from the thermal dilepton and Drell-Yan contributions. At RHIC e/i coincidence mea- 



surements could prove useful. Charm was first measured by this method at the ISR J26 
and such coincidence measurements are planned at RHIC [p7fl . Opposite sign e/i pairs 
from decay of two charmed particles are correlated if only one DD pair is produced in an 
event, which is the case for pp and light ion collisions at RHIC. If a = 1 in eq. (4), up 
to five primary DD pairs could be produced in Au+Au collisions, yielding uncorrelated 
e/i pairs. The production rate of Drell-Yan and thermal dilepton pairs is much less than 
one pair per event so uncorrelated e/i production from these processes is very unlikely. 
However, thermal charm may be important if the thermalization time is short, and cause 
significant additional production of uncorrelated e/i pairs. Pairs of like-sign electrons may 
offer a better measure of charm production if more than two DD pairs are produced per 
event. This is probable for all ion species in the LHC environment, where e/i coincidence 
measurements will not be possible in any case ^ffl. 

Figs. 5 and 6 suggest a very large correction to the measured dilepton distributions if 
subtraction of the charm contribution is attempted. However, these figures do not include 
the experimental acceptance, smaller for leptons from charm decay than for Drell-Yan or 
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thermal dileptons. Table 1 shows the acceptances in the muon channel of the RHIC 
PHENIX detector at M = 2, 4, and 6 GeV and y = 2. The Drell-Yan and thermal 



acceptances are from four to twenty times larger than for the charm decay leptons [27 
Consequently, the thermal and Drell-Yan contributions do not represent only a few percent 
of the total dilepton yield and subtraction of the charm contribution may thus be feasible. 
This is illustrated in Table 2 which shows the ratios of the cross sections without, then 
with the effects of acceptance. It is important to note that in the case of dileptons from 
uncorrelated charm decays, the acceptance drops by another factor of nearly two due to the 
larger rapidity difference between the leptons. It should be possible to optimize a rapidity 
cut to further discriminate against dileptons from uncorrelated charm. Whether such a 
cut can be used to reject a majority of the combinatorial lepton pairs from uncorrelated 
charm can only be addressed for specific detector geometries. 

If it is possible to reliably remove the charm decay contributions to the dileptons, 
the significant differences between the rapidity distributions of thermal and Drell-Yan 
dileptons could be measured. The Drell-Yan contribution grows slightly with y while 
the thermals decrease with y. The correlation of thermal production with the hadronic 
multiplicity and the effect of shadowing on Drell-Yan should lead to opposite trends in 
the mass number dependence of these two contributions. Stronger shadowing for higher 
mass nuclei favors the larger y region for Drell-Yan pairs. Increasing the mass number 
should narrow the rapidity distribution due to the stronger stopping and consequently 
also narrow the thermal dilepton distribution. The shape of the distributions in the 
central several units of rapidity should thus provide an indication of the dilepton source. 
A significantly broad rapidity and mass coverage at RHIC and LHC would allow such a 
study. 

In summary, we have calculated the rapidity distributions of dilepton production from 
various sources in heavy ion collisions at RHIC and LHC. The calculated distributions 
cover four or more units of rapidity, suggesting the contributions to dilepton measure- 
ments over a wide range of rapidities. A large experimental rapidity coverage is desirable 
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to allow a search for variations in the shape of the dilepton rapidity distributions in pp, 
pA and A A collisions. Charmed particle decays into leptons are very important, but the 
kinematics of the charm decays differ from that of thermal or Drell-Yan dilepton pro- 
duction. Consequently the acceptance of realistic detectors makes the charm background 
which must be subtracted to study the thermal dileptons more tractable. Furthermore 
the significant production of thermal charm may offer a novel way to study the collisions, 
since the rate is comparable to direct thermal dileptons. Experiments could search for 
broadening of the dilepton distribution in AA collisions after fixing the level of charm 
production in pp interactions. 

We would like to thank S. Gavin and G. Young for useful discussions. RV thanks Los 
Alamos National Laboratory and the University of Jyvaskyla for hospitality. 
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Mass 


2 GeV 


4 GeV 


6 GeV 


Drell-Yan 


3.9% 


5.4% 


4.6% 


Thermal dileptons 


4.9 


6.8 


5.8 


charm decays 


0.88 


0.89 


0.25 



Table 1: Acceptance in PHENIX muon arm for dileptons from various sources at y — 2. 



Ratio of dilepton source to Drell-Yan 



Mass 


2 GeV 


4 GeV 


6 GeV 


Thermal (C=l) 


6.1 


0.25 


0.013 


Charm 


66. 


33. 


3.8 


Thermal Charm 


14. 


0.73 


0. 



Acceptance correction included 



Mass 


2 GeV 


4 GeV 


6 GeV 


Thermal (C=l) 


7.7 


0.31 


0.016 


Charm 


15. 


5.5 


0.21 


Thermal Charm 


3.1 


0.12 


0. 



Table 2: The upper section gives the ratio of dileptons produced by other sources to 
Drell-Yan production at RHIC. The lower section shows the effect of the PHENIX muon 
acceptance on these ratios at y — 2. 

Tables 



15 



Figure Captions 



Figure 1. Rapidity distributions of Drell-Yan pairs at M = 2 (solid), 4 (dashed), and 6 
(dot dashed) GeV in a central Au+Au collision at RHIC (a) and the LHC (b). 

Figure 2. Initially produced charmed pair distributions (solid curves) from a central 
Au+Au collision at RHIC (a) and at the LHC (b). Also shown are the lepton pair dis- 
tributions resulting from semileptonic decays of charm at M D -^ = 2 (dashed), 4 (dot 
dashed), and 6 (dotted) GeV. The lepton pairs from the decays have been scaled up by 
a factor of 100. Thermal charm production is shown in (c) at RHIC and (d) at the LHC. 
The decay leptons from thermal charm at RHIC are scaled up by a factor of 1000 for 
M = 2 GeV (dashed), 10 4 at M = 4 GeV (dot dashed), and 10 5 at M = 6 GeV (dotted). 
At the LHC, the decays with M = 2 GeV are scaled up by 100 while those with M = 4 
(dot dashed) and M = 6 GeV (dotted) are scaled up by 1000 and 10000. 

Figure 3. The rapidity distribution of thermal dileptons produced at RHIC with a — 3 
in eq. (5) are shown for (a) and (b) M — 2 GeV, (c) and (d) M — 4 GeV, and (e) and 
(f) M = 6 GeV. The thermal curves represent the three choices of and T«, Tj = 1 fm 
(solid), C = 3 (dashed), and C — 1 (dot dashed). The full evolution (plasma, mixed 
phase, and hadron gas) is given in (a), (c), and (e) while the plasma contribution alone 
is shown in (b), (d), and (f). 

Figure 4. The same as in Fig. 3 for the LHC with a — 5. 

Figure 5. We compare the contributions to the dilepton spectrum at RHIC for (a) M = 2 
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GeV, (b) M = 4 GeV, and (c) M = 6 GeV. Drell-Yan (dashed), direct thermal dileptons 
with C = 1 (solid), thermal (dotted) and initial (dot dashed) charm production are in- 
cluded. 

Figure 6. The same as in Fig. 5 for the LHC. 



17 



